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' I. INTRODUCTION 

Measurements of time-dependent CP asymmetries in B° decays provide information about the 
irreducible phase contained in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1], 
which describes CP violation in the Standard Model. If a specific B decay final state has con- 
tributions from more than one amplitude and these amplitudes have different CP- violating weak 
phases, interference can produce a non-zero CP asymmetry. An essential ingredient in CP viola- 
tion measurements in B decays is flavor tagging. In this paper, we point out a subtlety of flavor 
tagging that has been overlooked or ignored in most recent CP violation analyses, describe the 
impact of this omission, and propose how to address it in some future measurements. 

In the current asymmetric i?-factories [2], PEP-II and KEKB, B°B° meson pairs are produced 
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in e + e~ interactions at the T(4S) resonance, where the pair evolves coherently in a P-wave state 
until one of the B mesons decays. Typically, one B decay is fully reconstructed and the flavor 
(whether it's a B° or B°) of this B, at the time of the other P's decay, is inferred from the decay 
products of the other B (the tag B). At the time of the tag B meson decay, the B mesons are 
known to be in opposite flavor states. In terms of the time difference between the two B decays, 
At = i rec — t tag , the time-dependent CP asymmetry is defined as 



where N is the number of events at At with a B° or B° as the tag B. 

Charged leptons and kaons are often used to infer the flavor of the tag B meson. The charge of 
a lepton from a semi-leptonic B decay has the same sign as the charge of the b quark that produced 
it. For example, a high-momentum e + (e~) would indicate that the tag B was a B° (B°) at the 
time of its decay. Similarly, a K + (K~) more often than not comes from a B° (B ). This works 
because the most likely b decay is b — > c and the most likely c decay is c — > s; thus the s quark 
usually has the same charge as the b quark. The lepton or kaon charge does not always correctly 
indicate the tag-P flavor. Mistags can come from incorrect particle identification or other B decay 
chains that produce wrong-sign leptons or kaons. The mistag fraction must be measured in order 
to determine the true CP asymmetry from the measured one. 

It is usually assumed that the measured CP asymmetry is entirely due to the interfering ampli- 
tudes contributing to the fully reconstructed B decay mode, and that the individual tagging states, 
such as P° — > D + ir~, are dominated by a single B decay amplitude. In other words, if only one 
B decay amplitude contributes to the tagging final state, it is safe to assume that all interference 
effects, such as CP violation, are due to the evolution of the fully reconstructed B. This assump- 
tion, which is valid for semi-leptonic B decays, ignores the possibility of suppressed contributions 
to the tag-side final state with different weak phases, such as happens for non-leptonic decays. 

These suppressed contributions may be important for kaon tags. For example, the D + ir~ final 
state with D + — ► K~it + it + , which is usually associated with a P° decay, can also be reached from 
a P° through a b — > ucd decay. Its amplitude is suppressed relative to the dominant P° decay 
amplitude (6 — > cud) by a factor of roughly \{V* b V cd ) / {V^ V* d )\ ~ 0.02, and has a relative weak 
phase difference of 7. Both Feynman diagrams are shown in Fig. 1. The tag-side b — > cud and 
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b — ► ucc/ amplitudes interfere, and, through the coherent evolution of the B°B° pair, alter the time 
evolution of Acp(At). The subject of this paper is to investigate the consequences of this small 
tag-side interference in some of the standard time-dependent CP-asymmetry measurements at B 
factories that use coherent B decays. 




CKM-favored amplitude doubly-CKM suppressed amplitude 



FIG. 1: The CKM-favored amplitude (left) and doubly-CKM-suppressed amplitude (right) for the final 
state D + iv~ . With respect to the dominant contribution, the latter is suppressed by the approximate ratio 
\(VubVcd) / (V c bV* d )\ w 0.02 and has a relative weak phase difference of 7. 

In Sections II - VI, we review the general formalism for describing the coherent evolution of the 
B°B° system, define our notation for describing the tag-side amplitude, and state the assumptions 
we employ in our analysis. In Section VII A, we evaluate how tag-side interference affects the 
mistag fraction measured from the amplitude of the time-dependent mixing (not CP) asymmetry. 
We find that the tag-side interference effects are not simply absorbed into the mistag fractions and 
that, to first order, the mistag fractions are unchanged by tag-side interference. In Section VII B, 
we evaluate the uncertainty, due to tag-side interference, in the standard mixing-induced CP 
asymmetry measurements - sin 2(3 from J/ij;K s and the CP asymmetry in 7r + 7r~. We find that 
the uncertainties are at most 5%, in the most conservative estimation, and can be limited to < 2% 
in most cases with reasonable assumptions. Finally, in Section VIII, we evaluate how tag-side 
interference affects some of the time-dependent techniques that have been proposed for measuring 
7 (e.g. the time-dependent analysis of D* + tt~). Here, we find that tag-side interference effects 
can be as large as the signal asymmetry. We propose a technique for performing the analysis in a 
general way, which does not require assumptions about the size of tag-side interference effects and 
maximizes the statistical sensitivity to (2(3 + 7). We summarize our conclusions in Section IX. 
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II. GENERAL COHERENT FORMALISM 

In this section, we define our formalism for describing the time evolution of a pair of neutral B 
mesons that are coherently produced in an Y(AS) decay and then subsequently decay to arbitrary 
final states ft and / r at times tt and t r , respectively, measured in the parent B meson's rest frame. 
The "t" ( "r" ) subscript refers to the tag (reconstructed) B meson or its final state. The amplitude 
for this process is proportional to 

A = (ft\B phy s(tt))(fr\B° phys (tr)) ~ (ft\Bl hy M))(fr\B° phy s(tr)) , (2) 

where B^ hys (t) (B phys (t)) denotes an initially-pure B° (B°) state after a time t. The relative minus 
sign between the terms reflects the antisymmetry of the P-wave B°B° state. Integrating over all 
directions for either B and the experimentally-unobservable average decay time (tt + t r )/2, we 
obtain a corresponding decay rate proportional to (At = t r — t t ) 

F(At) = e- r l A *l \a + g + (At) + a_ 5 _(At)| 2 , (3) 

where T is the average neutral B eigenstate decay rate and we define 

g±(At) = - (V iAmA */2 e -ArAt/4 _j_ e +iAmAt/2 e +ArAi/4^ ^ 

in terms of the differences between the eigenstate masses (Am) and decay rates (Ar). 

The time-independent complex parameters a± in Equation (3) can be written generally as 



a + = A t A r — At A r , a_ 



-Vl-z 2 (J A Ar - V -At a) + Z (At Ar + At Ar) , (5) 



where Ak (Ak) is the B° (B°) decay amplitude to fk- The complex ratio q/p parameterizes 
possible CP and T violation (\q/p\ ^ 1) in the time evolution of a neutral B state, while z, which 
is also complex, parametrizes possible CPT and CP violation (z ^ 0) in the time evolution. Note 
that exchanging the r and t subscripts changes the overall sign of a+, <7_, and At, leaving Eq.(3) 
unchanged, which is required since the distinction between the B that is reconstructed and the B 
that is used for flavor tagging is arbitrary at this point. Explicitly, we are using the conventions 



q = _ / M* 2 -TV2 

p V M12 - iT 12 /2 ' 1 ' 



where M and T are the hermitian matrices of the effective Hamiltonian. The eigenstates of the 
effective Hamiltonian are defined as 

\B L ) = p\B )+q\B ) (7) 

\B H ) = p \B°) - q \B°) , (8) 

and Am = ran — rriL, which is positive by definition. If z = 0, as expected in the Standard Model, 
the two terms a±g±(At) in Eq.(3) describe the cases where the surviving meson undergoes a net 
oscillation (— ) B° <-> B° or not (+) between t t and t T . Combining Equations (3-5) we obtain 

F(At) = e" r| Atl [R cosh(ArAt/2) + C cos(AmAt) + S' sinh(ArAt/2) + S sin(AmAt)] (9) 

with coefficients which satisfy the constraint C 2 + S 2 = R 2 — S' 2 , and are given by 

R= - (|a+| 2 + |a_| 2 ) , S'= -Re(a* + a-) , 

C= i (|a+| 2 - |a_| 2 ) , S = +lm(a* + a^) . (10) 

In the following, we assume CPT invariance so that z = 0, and moreover we take AT/T <C 1. 
Thus the term S' no longer enters and cosh(ArAt/2) is replaced by unity. The Standard Model 
predicts \q/p\ — 1 = (2.5 — 6.5) x 1CT 4 [3], so we will assume \q/p\ ~ 1. The resulting time 
dependence, when the tagged meson is a B°, is 

F(At) = e" r l A *l [R + C cos(AmAt) + S sin(AmAt)] , (11) 

and correspondingly when the tagged meson is a B° 

F(At) = e~ r|A ' 1 [R + Ccos(AmAt) + Ssm(AmAt)] . (12) 

III. CHARACTERIZATION OF THE TAGGING AMPLITUDE 



The strength of the doubly-CKM-suppressed (DCS) decays can be expressed in terms of the 
traditional parameter [4] 
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This combination is independent of the choice of phases for the B° and B° states. Suppose |/) is 
a final state that is ostensibly the result of a B° decay. For example, if |/} represents the tag B, 
a K + would indicate that the tag B decayed as a 5°, assuming the dominant b — > cud transition 
occurred. Then 

\ f = rf e- 2i ^e iS f , (14) 

where 77 is a real number of order 0.02 and 5f is the strong phase difference of the B° decay relative 
to that of the B° decay, assuming b — > cud and b — > ucd transitions for the B° and B° decays 
respectively. If, for this final state, there is only one mechanism contributing to the B° decay and 
to the B° decay, then for the CP conjugate state |/) we have 

X f = — e -w-i-y e - iS f . (15) 
r f 

We shall make the assumption of a single contributing amplitude except as noted below. 

Because the DCS amplitudes are only about 2% of the allowed amplitudes, in what follows we 
shall drop all terms that are quadratic or higher in this suppression. In practice we combine many 
final states / in a single tagging category, / G T. For the tagging category we then have effective 
values of r' and 5' defined by 

E / £ t tf\A f \ 2 r f e iS < 
where e/ is the relative tagging efficiency for the state /. Notice that 



so there is a tendency for contributions from different tagging states to cancel, unless all contri- 
butions have nearly the same strong phase. Equation 16 holds only if terms of order r 2 can be 
ignored, as we are assuming. 



IV. TIME-DEPENDENT ASYMMETRY COEFFICIENTS 

In this section, we evaluate the coefficients R(R), C(C), and S(S) of Eqns. 11(12). There are two 
specific cases that we will consider - the "mixing" case, where the reconstructed B meson decays 
in an apparent flavor eigenstate (e.g. D* + ir~, normally assumed to originate from B° decay), 
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and the "CP" case, where the reconstructed B has decayed into a CP eigenstate. Dropping a 
common factor At A T (p/q), we can write a + and a_ in terms of the A parameters for the tag and 
reconstructed B mesons as 

a + = \ t — A r 

a_ = l-A t A r . (18) 

Quite generally then, 

Kl 2 = |A t | 2 - 2 Re A t A* + |A r | 2 
\a_\ 2 = 1-2 Re A t A r + |A t | 2 |A r | 2 
Imo; fl _ = ImA r (1 - |A t | 2 ) - Im A t (1 - |A r | 2 ) . (19) 

Table I gives the coefficients for the mixing case, where for the reconstructed B meson final state 
we have dropped the subscript / from the amplitude ratio r and from the strong phase difference 
5 in A r , defined by Eq.(14). The only deviation from the familiar case with no DCS contributions, 
to first order in r and r', is the presence of a small S(S) coefficient. Figure 2 shows an illustration 
of the time evolution for when the flavor of the two B mesons at the time of decay was opposite 
(unmixed) or the same (mixed). The nominal (r = r' = 0) case is contrasted with an example of 
a non-zero DCS contribution in the reconstructed B amplitude and with an example of non-zero 
DCS contributions to both the tag and reconstructed B amplitudes. The amplitude ratios r and r' 
have been enlarged by x5 with respect to the expected value (0.02) so that the DCS contributions 
are more clear. 

Table II gives the coefficients for the CP case. All three coefficients receive corrections linear in 
r' . Figure 3 is an illustration of the the corrections to the time evolution for B° and B° tagged CP 
events, also with the DCS amplitude ratio r' enlarged by x5 to make the differences more visible. 

V. COMPLETELY INCLUSIVE TAGGING CATEGORIES 

We can relate the effective r' and 5' to the 2x2 matrix T that generalizes the decay rate for 
the B°B° system. Let T be the class of states DX, where X represents non-charmed hadrons. 
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[tag=B°(A"+), rec=B° ] [tag=B°(K+), rcc=B° ] [tag=B°(A _ -) 1 rec=£° ] [tag=B°(A"-), rec=B° ] 

^ t r l e -2ifj-i-,+iS' r l e -2i0-i~/+i6' l re -2i0-ii-iS' J Te -2if3-i~/-i6' 

^ re — 2i/3 — i^+iS l e — 2ifl — i~f — i5 re — 2ij3 — i-y + i5 1 g — 2i/3 — 27 — iS 

r r r 

k+l 2 4r ^ 

l«-| 2 1 

-rsin(2/3 + 7 -5) sin(2/3 + 7 - 6') -fc sin(2/3 + 7 - 5) - ^7 sin(2/3 + 7 + 

+r'sin(2/3 + 7-5') -± sin(2/3 + 7 + 5) + ^8111(2/3 + 7 + 5') +^ sin(2/3 + 7 + 5) 

_R 1 1 1 1 

C -1 1 1 -1 

5 -2rsin(2/3 + 7-5) -2r' sin(2/3 + 7 - 5') 2r sin(2/3 + 7 - 5) -2r' sin(2/3 + 7 + 8') 

+2r' sin(2/3 + 7 - 5') -2r sin(2/3 + 7 + 5) +2r' sin(2/3 + 7 + 5') +2r sin(2/3 + 7 + 8) 



TABLE I: Contributions to the time dependence of tagged decays when the reconstructed decay is an appar- 
ent flavor eigenstate, with doubly-CKM-suppressed decays considered only to first order. The dependences 
proportional to 1 and to cos AmAt are unaffected. A small sin AmAt term is induced. Appropriate factors 
of r and r' have been removed to scale R to unity. 

Neglecting the relative tagging efficiency ej for the moment, we have 

J2 -A* f A f = \Af\ 2 \f = J2( B °\ n \f)(f\ n \ B °) r f e ~ 2lP ^ +t5f = T DX r'e~ 2l ^ +iS ' , (20) 

/GT P feT /GT 

where Tdx is, up to a trivial normalization, the partial width of B° into the class of states of the 
form DX. On the other hand, we can write 

E l A *i*f = E l^\n\f){f\n\W) = f r M12 , (21) 

TgT 1 FgT p p 

where Fnx 12 is the contribution of states of the form DX to the off-diagonal part of the V matrix. 
So 



I~2if3-ij+i8' _ Q r /-p 

re = -1 ox 12/I 

P 



(22) 
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[tag=5° (K+ ) , rcc= B c P ] [tag=B° (K " ) , rec=B CP ] 

^ t r l e -2i(3-i-y+i&' I_ e -2il3~ii-i8' 

A r A C p Acp 

\a+\ 2 \X CP \ 2 - 2 Re r^-atf-fry-H*' A^ p ^-2iRc A CP e 2 ^+ J T+ l5 ' 

|a_| 2 l-2Rcr'e- 2 ^-^+ 4 ' 5 'Acp - 2^ Rc \ C pe- 2 ^-^- t5 ' 

Ima; a _ Im A C p +r'(l- |A£ P |) sin(2/3 + 7 - <*') Im A CP + £(1 - |A CP | 2 ) sin(2/3 + 7 + £') 

i? 1+l f" pl - 2r' Rc A CP cos(2/3 + 7 - <*') - 2r' Re A C p cos(2/3 + 7 + 5') 

C I a cpM + 2r ' Im A CP sin(2/3 + 7 - 5') 1 ~'^ pl + 2r' Im A CP sin(2/3 + 7 + 5') 

5 ImA C p + r'(l- |A 2 7P |)sin(2/3 + 7 -5') - Im A C p + r'(l - |A£ P |) sin(2/3 + 7 + 5') 



TABLE II: Contributions to the time dependence of tagged decays when the reconstructed decay is a CP 
eigenstate, with doubly-CKM-suppressed decays considered only to first order. Appropriate factors of r and 
r' have been removed to scale R to unity in the limit in which the doubly-CKM-suppressed decays vanish. 



If tagging does not capture every state, we can think of T\2 and T as effective quantities, limited 
by the partial sum over states. However, if that sum were complete, then 5' would vanish. To see 
this, imagine using as a basis of states fs not the physical states that are observed but instead a 
basis of states that are eigenstates of the S matrix, that is, a basis of states that each scatter into 
themselves. Because we are summing over all states in a collection connected by strong interactions, 
there is such a basis. Then the final state interaction phases associated with Af s and Af s would 
both be e lSf s . These would cancel in A^ s Af s . In general, because tagging is incomplete, we cannot 
assume that 5' vanishes. In reality, the relative tagging efficiency e/, set to one in Eq.(20), is not 
the same for all of the states in T, so the tagging category representing the class of states DX is 
not completely inclusive. 
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FIG. 2: Time-dependent decay distributions for the final state D*~-k + , for a) a B° tag, and b) a B° tag. 
(2/3 + 7) is set to the value 1.86. The situation with no doubly-CKM-suppressed contribution on both the 
tag-side and reconstruction-side is indicated with the solid line. The dotted line has r = 0.1 and 5 = 0, but 
no tag-side interference. The dashed line represents the example with r = r' = 0.1, 6 = 0, and 5' — ir. In 
these examples, the r and r' values are x5 the expected values in order to clearly illustrate the differences 
with respect to the case with r = r' = 0. 



VI. ESTIMATED SIZE OF DOUBLY-CKM SUPPRESSED AMPLITUDE 

In the Introduction, we gave an estimate for the size of the DCS amplitude (r), relative to 
the favored amplitude, to be approximately 0.02, which is simply the ratio of the CKM elements 
involved, \(V* b V cd ) / (V cb V* d )\. Here, we discuss the uncertainty of this estimate as well as what 
can be assumed, if anything, about the strong phase difference (5) between the DCS and favored 
amplitudes. 

We use measured charm branching fractions as a test of our simple amplitude ratio estimate. The 
charm decay D° — > K + ir~ is doubly-CKM suppressed relative to the favored D° — > K~tt + decay. 
The amplitude ratio prediction, based solely on the CKM elements, gives r ~ |(^d^us)/(^cs^ud)l ~ 
0.048 . The experimental value from the branching fractions [5] is 0.062±0.005, which is within 25% 
of the ratio of CKM elements. For the singly-CKM suppressed charm decays D° — > K + K~ and 
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— > 7r + 7r we would estimate amplitude ratios relative to the allowed amplitude of IV^/V^I « 
IWKsl ~ 0.23, while the branching ratios give 0.329 ± 0.007 and 0.194 ± 0.005 for K + K~ and 
7r + 7r~, respectively. 

The decay B° — > D + ir~ is doubly-CKM suppressed, but this branching fraction has not been 
measured. We can estimate its branching fraction from the related decay mode B° — > Dfir~, 
which has been observed recently [6], and has a branching fraction of (3.2 ± 0.9 ± 1.0) x 10~ 5 . The 

1.0 
0.5 

0.0 

1.0 

0.5 

0.0 
0.5 
0.0 
-0.5 

-10 -5 5 10 

At (ps) 

FIG. 3: Time-dependent decay distributions for the CP eigenstate J/tpKg, with a) a B° tag, and b) a B° 
tag. We set (2/3 + 7) to 1.86. The situation with no tag-side interference is indicated with the solid line. 
The dotted line represents the case with r' = 0.1 and 5' = 0, and the dashed line has r' = 0.1, and 5' = n. It 
should be noted that, adding a non-zero DCS contribution, the slope and amplitude of the time-dependent 
asymmetry work in opposite directions. In these examples, the r' value is x5 the expected value in order 
to clearly illustrate the differences with respect to the case with r' = 0. 
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amplitude ratio for B° — ► D + 7r , relative to B° — > D 7r+ is estimated to be 





where we have used /d//d s = 1-H ± 0.01 ± 0.01 (from [7]) to approximate S77(3) breaking effects. 
This is in good agreement with the naive estimate of 0.020, albeit with a large uncertainty. 

There are some theoretical arguments for expecting the strong phase difference 5 to be small [8], 
but we know of at least one case where a non-trivial strong phase has been observed in B decay. 
The strong phase difference between the longitudinal and parallel polarization amplitudes of the 
transversity basis in B — > J/ipK*(892) has been measured [9] to be 2.50 ± 0.22, which is about 
3a from it, in contradiction with the factorization prediction of or it. The size of the effective 
amplitude ratio (r'), given by Equation 16, depends on the 5 values of the final states included in 
the tagging category. As Eq. 17 shows, varying 5 values between the states will tend to reduce r' . 

Given the « 50% uncertainty on the DCS amplitude ratio r for individual final states and 
the general lack of knowledge concerning strong phase differences, we conclude that the most 
conservative assumptions regarding the effective parameters r' and 5' would be to allow r' values 
from (full cancellation in the sum) up to 0.04 (no cancellation with some enhancement over our 
0.02 estimate) and to allow any value of §'. 



In this section, we will discuss the uncertainties due to tag-side interference on some common 
time-dependent asymmetries. In addition to the assumptions that we have already made (i.e. 
z = 0, Ar/r = 0, and \q/p\ = 1), one usually assumes that the tag-side amplitude is dominated 
by a single contribution, or r' = 0. The time dependent coefficients in Eqns. 11 and 12 simplify 
considerably with this assumption. For the case where the reconstructed B is a CP eigenstate, we 
have 



which can be seen from Table II with r' set to zero. For the case where the reconstructed B is in 
an apparent flavor eigenstate, the coefficients in Table I with r' = give 



VII. UNCERTAINTIES IN FITTED ASYMMETRIES 



Rcp = Rcp 



Ccp = —Ccp 



Scp = -Sep , 



(23) 



R 



mix 



unmix 



Cmix — Ci 



unmix 




(24) 
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where the "mix" ( "unmix" ) subscript refers to the case where the tag and reconstructed B mesons 
were the same (opposite) flavor at the time of decay. In the rest of this Section, we will evaluate 
the bias on the fitted coefficients when fitting the data with the assumptions in Eqns. 23 or 24 of 
nonzero tag-side interference. 

In the relations above, the R coefficients are independent of the final state configuration, so 
they are usually absorbed into the C and S coefficients by fitting for C = (C/R) and S = (S/R). 
A fairly reliable estimate of the fitted C coefficient is simply the asymmetry at At = 0. This would 
be 



for a CP asymmetry, or 



C+R-C-R 
Cfit ~ C+R+C+R (25) 



n Cunmix -^unmix Cmix Rmix /„„\ 

Cfit ~ r Ir Tr — Ir — ^ lb > 

'-'unmix ~r "unmix < '-'mix ~r -fornix 



for a mixing asymmetry. A similar, but slightly less reliable, estimate for the fitted <S coefficient 
in a CP asymmetry is simply the flavor-averaged S coefficient, or 

* «\ (| - |) • («) 

Precise estimates can be derived using a simple maximum likelihood technique, where the 
likelihood to be maximized with respect to C nt and >S nt is 

/oo 
dAt e- r|A *! [F(At) In F fit (At) + ~F(At) lnF fit (At)] , (28) 
-oo 

with Ffit and F& t evaluated using the assumptions in Eq. 23. We confirmed that Eqns. 25 and 27 
give reasonable estimates of C nt and <S nt with unbinned maximum likelihood fits of simulated data 
samples. 



A. Mistag calibration with flavor oscillation amplitude 

As was mentioned above, the sign of the tagging kaon charge does not always give the correct 
flavor tag. For example, CKM-suppressed D decays, such as D + — > K + K°, can produce wrong-sign 
kaons. Pions, incorrectly identified as kaons, can also produce wrong-sign kaons. The amplitude 
of any measured asymmetry using kaon tags will be reduced by a factor of (1 — 2uj), sometimes 
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called the dilution factor, where u> is the fraction of tagging kaons that have the wrong sign (mistag 
fraction). The mistag fraction u is usually measured from the amplitude of time-dependent flavor 
oscillations in a sample of reconstructed B° decays to flavor-specific final states [10]. The measured 
value of C will be a direct measurement of (1 — 2u>), which can then be used to translate measured 
CP asymmetry coefficients. 

To first order in r and r', the R and C coefficients are the expected ones, as can be seen in 
Table I. The only effect is in the S coefficient, which is usually assumed to be zero in the analysis of 
mixing data. This means that the measured mistag fractions will be unaffected by DCS amplitude 
contributions, either on the tag side or the reconstructed side, since our estimator for Cg t only 
depends on the R and C coefficients. Contrary to what one may guess, the corrections due to DCS 
amplitude contributions are not simply absorbed into the mistag fractions. 

Using Monte Carlo pseudo-experiments, we also find that Am^ is unaffected to the level of 0.001 
ps _1 if allowed to float in the fit. 



B. Fully reconstructed CP eigenstates 

The size of CP asymmetries in B decays to CP eigenstates are in general of order one in the 
Standard Model. For example, CP asymmetry in B — > J/ipKg (and related charmonium modes) 
has been measured to be Sa t = 0.735 ± 0.055 [11, 12]. Any deviations due to tag-side interference 

0.02) will be comparatively small (see Fig. 3), and can be treated as perturbations on the usual 
measurements. 

In what follows, the nominal values for the fitted CP asymmetry coefficients without any tag- 
side interference from doubly-CKM suppressed decays are defined as 

\X C p\ 2 - 1 , s 

Co ~ Jx^FTT (29) 

So = r^fj- (30) 

The expected fitted coefficients, when the fit is performed with the assumptions in Eq. 23, can be 
found by inserting the R, C, and S values from Table II into Eqns. 25 and 27. Working to first 
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order in r', we find 

Cfi t = C [1 + 2r' cos 5' {Q cos(2/3 + 7) - S sin(2/3 + 7)}] 

- 2r sin5 / {5 cos(2/? + 7) + £ sin(2/3 + 7)} (31) 

<S fit = So [1 + 2r cos 5'G cos(2/3 + 7)] + 2r sin <5'C cos(2/3 + 7) , (32) 

where Q = 2ReAc*p/(|Ac*p| 2 + 1). Note that, with respect to the nominal values, there are both 
multiplicative and additive corrections which are proportional to cos 5' and sin 5' respectively. In 
the limit of a vanishing effective tag-side strong phase difference (8' — > 0), only the multiplicative 
corrections remain. 

For B° — > J/ipKg, the dominant tree and penguin amplitude contributions share the same 
weak phase. The highly suppressed u-quark penguin, which has a different relative weak phase, is 
typically ignored, giving the Standard Model prediction of Aj/^o = — e _l2/3 . Inserting this into 
Eqns. 31 and 32 gives 

C &t [J/ipK^] = -2/ sin 7 sin y (33) 
S m [J/ijK s ] = S [1 -2r'cos(5 / {cos2/3cos(2/3 + 7) + /Csin2/?sin(2/3 + 7)}] , (34) 

with Cq = and So = sin 2/3. The last term in Eq. 34 proportional to K. is a correction to the 
simple estimate given by Eq.(32). The correction /C was derived from the more precise likelihood 
analysis given by Eq.(28). The value of K, is between 0.10 and 0.35, depending on the value of 
sin 2/3. If we assume sin 2/3 = 0.74 and allow 7 to be in the range [45°, 90°], then JC = 0.28 and 
the magnitude of the deviation of <Sfit away from the nominal value So is < 0.7 r. The size of the 
deviation of CstlJ/ipKg] could be as large as 2 r' . These corrections to <Sfit = So and Cat = could 
be as large or larger than Standard Model corrections [13]. 

The uncertainty estimates in the previous paragraph apply to a measurement that only uses 
kaon tags. In practice, all useful sources of flavor information from the tag side B are employed 
in order to maximize the sensitivity of the measurement. The statistical error on the measured 
asymmetry scales as Qi, where each flavor tagging category contributes Qi = ej(l — 2uoi) 2 

and €i is efficiency for category i. Lepton flavor tags do not have the problem of a suppressed 
amplitude contribution with a different weak phase, so we assume that r' = for lepton tags. If 
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a measurement uses both lepton and non-lepton tags, the magnitude of the tag-side interference 
uncertainty will be scaled down by a factor of Qnon— lcp /(Qicp + Q non— lcp)- For example, the BaBar 
flavor tagging algorithm[ll] has roughly Q\ cp « 0.1 and Q non -lcp ~ 0.2. This gives a reduction of 
the tag-side interference uncertainty of about a factor of 2/3. 

The CP asymmetry for B — > 7r + 7r~ is more complex. This decay has both tree and penguin 
amplitude contributions which are comparable in magnitude, have different weak phases, and have 
an experimentally unknown relative strong phase difference. Equations 31 and 32 do not become 
more transparent after inserting the value for X W7r given below 



where the t-quark penguin has been absorbed into the tree and penguin amplitudes using unitarity 
of the CKM matrix, as in [14]. Clearly, both the reconstructed and tag B amplitudes now depend 
on 7, so care must be taken in evaluating the tag-side interference uncertainty, which in general 
can be as large as 2 r' for either the multiplicative or additive terms in Eqns. 31 and 32. 



One technique for measuring or constraining 7 is to perform a time-dependent analysis of a 
decay mode that is known to have a non-zero DCS contribution, such as D* + ir~ [15]. The time- 
dependent asymmetry coefficients are those given in Table I. In the usual case, tag-side interference 
is ignored {r 1 = 0) and the amplitude of the sin Am At term is 2r sin(2/3+7±(5), where r is the ratio 
of the DCS to CKM-favored amplitude contributions for the reconstructed, or non-flavor-tag, B 
and 8 is the strong phase difference between the two amplitudes. Measuring r and sin(2/3 + 7 ± 5) 
simultaneously is very challenging, so it is likely that r will have to be constrained from other 
measurements [6]. 

Since both r and r' are expected to be of the same order (« 0.02), it is clear that tag-side 
DCS interference can not be treated as a perturbation on the usual case. This effect is illustrated 
in Fig. 2. The time dependent analysis should be performed in a way that is general enough to 
accommodate r' ~ r and any value of 5'. 

Table III gives the sin Am At coefficients, taken from Table I, for the 4 combinations of recon- 
structed and flavor tag B final states, where we have neglected r 2 , rr', and r' 2 contributions. It is 




(35) 



VIII. MEASUREMENT OF SIN(2/3 + 7) WITH £>(*> 



7T 



17 



Symbol 




Reco 


Tag 


sin(AmAi) coefficient 


51 


B° 


(£)*-7T+) 


B° (#+) 


-2 r sin(2/3 + 7 - 5) + 2 r' sin (2/3 + 7 - 5') 


S2 


B° 


(£)— 7T+) 


B° (if-) 


2 r sin(2/3 + 7 - 5) + 2 r' sin(2/3 + 7 + 5') 


S3 


w 


(D* + 7T-) 


B° (K+) 


-2 r sin(2/3 + 7 + 5) - 2 r' sin(2/3 + 7 - S') 


S4 


B a 


(^* + 7T-) 


B° (K-) 


2 r sin(2/3 + 7 + 5) - 2 r' sin(2/3 + 7 + 6') 



TABLE III: The 4 coefficients of the sinAmAi term in the time-dependence of D*ir. The 2nd and 3rd 
columns give the interpretation of the observed final state (given in parentheses) in terms of the dominant 
amplitude. 



useful to rewrite the relations for the S coefficients in the following way 



51 = -a + b + c (36) 

52 = +a + b-c (37) 

53 = -a-b-c (38) 
SA = +a-b + c, (39) 

where the 3 variables to be determined in the time-dependent analysis are 

a = 2 rsin(2/3 + 7) cos 5 (40) 
b = 2 r' sin(2/3 + 7) cos 5' (41) 
c = 2cos(2/? + 7) (rshi(5 - r'sin<5') . (42) 

This parameterization makes no assumptions about the magnitude of r' or 5', and is attractive 
for several reasons. First, a does not depend at all on the tag-side parameters r' and 5'. In the 
case where 5 = 0, which is favored by some [8], a is exactly what one wants to know (sin(2/3 + 
7)). Secondly, this parameterization cleanly separates the flavor-tag symmetric and antisymmetric 
components; the a and c coefficients are diluted by a factor of (1 — 2uj), while the b coefficient is 
not, since it has the same sign for tag-side B° and tag-side B° events. The minimum number of 
independent parameters in which the S coefficients can be written is three. We recommend using 
the a, b, and c coefficients as the experimental parameters to be determined in the time-dependent 
asymmetry analysis. 
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The set of kaon tagging final states that yields correct tags is in general quite different from 
the set of final states that yields incorrect tags. This means that within a tagging category, the 
effective r' and 5' values for correct tags are different from those for incorrect tags. In the sum over 
correct and incorrect tags, the terms linear in r' that appear in the observables of the asymmetry 
are 

(1 - 2u)r'e iS ' = (1 - uj)r' c e i5 'c - ur'^'i . (43) 

This equation gives effective r' and 5' parameters in terms of the mistag fraction oo, effective 
parameters for correct tags (r' c and 8' c ) and incorrect tags (r[ and 5^). This implies that, in order 
to have a completely general parameterization in the data analysis, each tagging category (kaon, 
lepton, slow pion, etc.) must have different effective r' and 5' parameters, and thus different b and 
c parameters, due to the dependence on the mistag fraction u. One particular case that is relevant 
for a kaon tag category is when = 0. In this case r' = r' c (l —cj)/(1 — 2a;), which means that the 
effective r' is enhanced by a factor of (1 — cj)/(1 — 2u). 

The experimental knowledge of 5 depends on c, so even though the a parameter does not depend 
on r' and 5', one does not avoid uncertainties due to r' and 5' in the analysis. The best way to 
reduce this uncertainty is to take advantage of the fact that lepton tags are immune to the problem 
(r' = 0). If the fit is performed with an independent c coefficient for lepton tags, q ep combined 
with the a parameter measured by all flavor tagging categories will help resolve 5 and thus (2/3 + 7). 

If r' and 5' are not constrained from other measurements, one must allow for values of r' and 5' 
that are consistent with the measured values of b and c. Since it is possible to have a measured set 
of a, b, and c parameters that are consistent with r' = when r' ^ 0, one must always consider all 
r' values between and r' max consistent with b and c, where r' max is the largest allowed single-final- 
state value. This point is illustrated in Figure 4. The uncertainty on (2/3 + 7) due to r' and 5' is 
maximal when a is small. In this case, the sensitivity to (2/3 + 7) is mostly from the c coefficient 
and one must rely on flavor tag categories that are known to have r' = 0, such as lepton tags. 

Using Monte Carlo pseudo-experiments, we perform a simplified study of the impact of DCS 
tag-side interference on a system with only two tagging categories: one for unaffected lepton 
tags, and the other containing kaon tags. The significance ratio of both categories is set to 
Qiep/Qnon-lep = 0-6- All tests use the realistic value of 0.02 for r and r' . Each category shares the 
same a parameter. The lepton category constrains q ep , and the kaon category fits b and c. All fit 
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FIG. 4: An example of three amplitude configurations that all give the same set of S coefficients. For 
each set (a, b, or c) the S coefficients are consistent with no tag-side interference (r' — 0), while this is 
only true in the first case. Each configuration is represented by four diagrams showing the addition of the 
reconstructed and tag-side amplitude vectors in the complex plane. The observable coefficient (5*) is the 
imaginary part, represented by the vertical band on the right side of each diagram. The parameters for the 
three configurations are: r = 0.02; (2/3 + 7) = 2.10, 1.87, 1.73; r' = 0.00, 0.02, 0.03; 8 = 0.30, -0.53, -0.58; 
and 5' = NA, 1.57, 1.57. 

parameters are unbiased, and conform to Gaussian distributions. Compared to the situation with 
no DCS contribution, having one tagging category and identical errors for its two parameters a 
and c, the statistical error on a is unchanged, and that on c\ ep has increased by a ratio compatible 
with ((Qnon-lep + QiepVQiep) 1 / 2 = 1-6- The parameters a and b show a 20% correlation, while all 
other correlations are smaller than 1%. 

One experimental strategy for reducing the uncertainties due to r' and 5' would be to constrain 
them by performing a time-dependent analysis of a flavor-specific final state that has no DCS 
contribution (r = 0), such as D* + l~v. For such a final state, the undiluted b coefficient is the 
same as for D* + ir~ and c now has r = 0. This information can be used to recover the (2/3 + 7) 
sensitivity in the c coefficients in the signal sample that was lost due to the lack of knowledge of r' 
and 5' . Another option would be to include in the analysis events for which it was not possible to 
determine the flavor of the tag, so-called untagged events. From Equations 36 through 39, one can 
see that the untagged 5 coefficient for a reconstructed D*~ir + (D* + ir~) is equal to 51 + 52 = b 
(53 + 54 = —b), thus untagged events provide a further constraint on b. 
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The measured a, b, and c coefficients for the various tagging categories and sampies can be 
combined by forming a x 2 using the measured parameters and the inverted covariance matrix. 
This assumes that the measurement uncertainties on the a, b, and c parameters are Gaussian. A 
constraint on (2/3 + 7) can be derived from the \ 2 by scanning the \ 2 vs (2/3 + 7) where for each 
(2/3 + 7) value the \ 2 ls minimized with respect to the unknown parameters 6, 5', and r' . If there 
are no external constraints on r' and 5', such as from the analysis of D* + l~V suggested above, 
the b and c parameters from non-lepton tags do not provide much information, since r' must be 
varied from its minimum value compatible with b to its maximum possible value (for example, see 
Figure 4). The non-lepton-tag b and c parameters still must be included in the time-dependent fit, 
but they are not very useful in the x 2 analysis. 

Figure 5 shows an example of the \ 2 procedure for a hypothetical measurement where (2/3+7) = 
1.86, r = 0.02, and S = 0.9. The measured values of the a, b, and c coefficients were set to the 
correct values, so the \ 2 is zero at the correct and degenerate solutions. The two plots in Figure 5 
illustrate two cases: a) one with non-zero tag-side DCS interference, and b) one without tag-side 
interference. In addition to the curve which allows for any value of r', labeled 'envelope', additional 
curves with fixed values of r' are included. The statistical errors correspond to a measurement from 
D*tt in roughly 450 fb _1 of i?-factory data from one experiment including a constraint from D*lv. 

Three important conclusions can be drawn from Figure 5. First, comparing envelope curves for 
the r' = 0.02 case a) to the r' = case b), the measurements give nearly identical constraints on 
(2/3 + 7). This means that the uncertainty on r' and 5' does not affect the measurement. The only 
degradation with respect to the situation with no tag-side interference is that, when not including 
D* + l~V in the analysis, the non-lepton-tag c parapeters no longer contain useful information. 

The second conclusion is that if r' is non-zero, the constraint on (2/3 + 7) can be better than the 
case where r' is zero. If so, the b and c parameters in the D*lv sample will in general be non-zero, 
and one effectively adds a measurement of (2/3 + 7) from the tag-side B. This can be seen most 
clearly from the symmetry between the tag-side and reconstruction-side within the definitions of 
a, b, and c in Eqns. (40-42). An extreme, unrealistic example is given by the solid line in case a) 
of Figure 5. It shows what the constraint looks like if r' would equal 0.02, and if that information 
were known precisely and were included in the analysis. 

Thirdly, the result for D*tt alone, after varying r' to arbitrarily large values, is equivalent to the 
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FIG. 5: Scans of % 2 from measured a, b, and c coefficients as a function of (2/3 + 7) illustrating two cases: 
a) one with non-zero tag-side DCS interference, and b) one without tag-side interference. The solid curve 
in both a) and b) was made with the true value of r' . The dashed curve labeled 'envelope' in figure a) 
encloses from below all \ 2 curves made with r' ranging from zero to arbitrarily large values. The envelope 
curve in b) coincides with the % 2 curve made with the true r' value of zero. The dotted curves in both a) 
and b) illustrate x 2 curves with incorrect values of r'. The input values were (2/3 + 7) = 1.86, r = 0.02, 
and 8 = 0.9. The tag-side parameters for a) were r' — 0.02 and 8' = 0.3. The measured values of the a, b, 
and c coefficients were set to the correct values. The statistics of the hypothetical measurement correspond 
to roughly 450 fb _1 of B factory data from one experiment, including a constraint from D*lv. There is a 
discreet ambiguity that gives exactly the same curves after adding 7r to the horizontal axis (2/3 + 7). 

X 2 curve constructed from only a and q ep . In other words, when not including the D*lu sample in 
the analysis, the b and non-lepton-tag c parameters do not contribute to the sensitivity to (2/3 + 7). 
Again, however, these degrees of freedom must still be included in the data analysis. 

IX. CONCLUSIONS 

Interference effects between CKM-favored b — > cud and doubly-CKM-suppressed b — » ucd am- 
plitudes in final states used for flavor tagging in coherent B°B° pairs from T(4S) decays intro- 
duce deviations from the standard time evolution assumed in CP violation measurements at the 
asymmetric-energy B factories. To our knowledge, the uncertainty introduced by this interference 
has been neglected in most B factory CP violation measurements published to date, with the 
exception of [11]. The uncertainties introduced in the sin 2(3 measurement in (cc)K° decay modes 
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and the time dependent analysis of the ir + ir~ final state are at most of the order of 5% and can 
be limited to < 2% in most cases with reasonable assumptions. 

In proposed measurements of sin(2/3+7) which explicitly use interference between CKM-favored 
and doubly-CKM-suppressed amplitude contributions in the final state that is reconstructed, such 
as D*tt, tag-side interference effects can be as large as the interference effects one is trying to 
measure. In any such analysis, the data must be analyzed in a way that is general enough to 
allow for tag-side interference effects. We have proposed a general framework for dealing with 
tag-side interference effects in sin(2/3 + 7) measurements. It is possible to achieve an experimental 
sensitivity to (2/3 + 7) similar to the originally proposed measurements. 
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